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We report on high-precision measurements of the temperature dependence of the plane- 63 Cu NQR 
line frequency ^q(Cu2) and the linewidth in normal and superconducting le O and ls O exchanged 
YBa2Cu4C>8. Whereas z/q(Cu2) passes T c very smoothly without a discontinuity either in value nor 
in slope, the linewidth increases in the normal conducting phase down to T c and starts to decrease 
sharply in the superconducting phase to finally resume its high-temperature value of the normal 
phase. 

There is a well discernible oxygen isotope effect on the vq (Cu2) temperature dependencies. The 
temperature dependence of vq (Cu2) is described by an empirical expression consisting of two parts: 
one related to the thermal expansion of the lattice and the other due to charge redistribution during 
the formation of new electronic structures in the CUO2 planes. From the fit to the experimental data 
we determine for the conjectured formation of new electronic structures an energy scale A( ls O) = 
188.0(1.6) K and A( ls O) = 180.0(1.6) K. This results in a partial oxygen isotope effect coefficient 
olvq = 0.42(11) which is larger than both the spin-pseudogap coefficient apa = 0.061(8) and the T c 
coefficient a Tc = 0.056(12) |[. 

PACS numbers: 74.72.Bk, 76.60.Gv 
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I. INTRODUCTION 



In the past few years very strong experimental evi- 
dence has been found for static or dynamic charge in- 
homogeneities in strongly correlated electronic systems 
in particular in high-T c superconductors where an ar- 
ray of self organized one dimensional structures known 
as stripes can appear and with them the possibility of a 
pseudogap |2| . Static charge redistributions might occur 
at the formation of new electronic structures. The nu- 
clear quadrupolc resonance (NQR) frequency (vq) being 
proportional to the electric field gradient (EFG) can serve 
as a very sensitive monitor of any such charge redistribu- 
tion. The plane-copper (Cu2) nuclear quadrupole reso- 
nance frequency vq(Gu2) in YBa2Cu40s seems to signal 
such an event. Besides the effect of thermal lattice ex- 
pansion, the temperature dependence of ^q(Cu2) reveals 
an additional effect possibly connected to the charge re- 
distribution in the Cu02 planes induced by some kind 
of a new electronic structure. The temperature depen- 
dence of vq(Gu2) is very unconventional, particularly 
striking is the change of the course around 190 K where 
a broad minimum in the frequency appears jjj. In the 
same temperature region other anomalies in nuclear mag- 
netic resonance (NMR) and NQR quantities, for instance 
in Knight shifts, line widths, and relaxation times oc- 
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cur which signal an electronic crossover [§, §. To learn 
more about the CuC>2 planes especially about a possible 
charge redistribution during the pseudogap formation as 
well as to have very precise vq(Cxx2) values frequently 
demanded in other Cu2 NMR and NQR studies we de- 
cided to implement the Cu2 vq and linewidth data by 
additional measurements on the very same material on 
which we recently determined by Cu2 spin-lattice relax- 
ation the oxygen isotope effect of the spin-pseudogap M . 
Therefore we performed between 5 K and 350 K very ac- 
curate measurements of the temperature dependence of 
the Cu2 NQR line frequency i^q(Cu2) and the linewidth 
in normal and superconducting 16 O and 18 O exchanged 
YBa2Cu40g. It is our anticipation that accurate knowl- 
edge of the influence of different oxygen isotopes onto 
vq (Cu2) and its temperature dependence will add clues 
to the understanding of high-T c superconductors. 



II. EXPERIMENTAL DETAILS 

Sample preparation, oxygen exchange process and 
sample characterization are described in Ref. @. The 
18 O content in the 18 O sample is 88% as determined by 
the weight loss of 18 O material after back exchange of 
18 O with 16 O. Room temperature x-ray measurements 
show a small difference in the lattice parameters of the 
two oxygen isotope samples. The a, b, and c lattice pa- 
rameters are 3.8411(1) A, 3.8717(1) A, 27.2372(8) A for 
the 16 and 3.8408(1) A, 3.8718(1) A, 27.2366(8) A for 
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the 18 O samples. For the detection of the NQR signal 
we used a standard NQR pulse spectrometer employing 
spin-echo technique and echo recording in quadrature. 
The magnitude of the complex Fourier transform of the 
whole echo delivered the shape and the position of the 
line. The two oxygen exchanged samples were inserted 
into a probe head with two identical resonance circuits 
which allowed simple switching of the electronics from 
one sample to the other thus minimizing the effects of 
any possible slow drifts in the characteristics of the equip- 
ment and temperature. To guarantee a proper excitation 
and detection of the relatively broad line we increased the 
damping of the resonance circuits by additional 7 Q resis- 
tances. Further we used 6 mm diameter coils and highest 
power available (1 kW) such that a pulse of only 1.8 /is 
duration corresponded to a ir pulse. We also took care 
to tune the resonance circuits as well as to excite the line 
with pulses at a frequency very close to the middle of 
the line, the difference between the two never exceeding 
more than 10 kHz. During the whole experiment we fol- 
lowed a rigid measurement procedure to ensure an equal 
treatment of the two samples. 



III. ANALYSIS AND RESULTS 



dence at T c with the subsequent narrowing of the line in 
the superconducting state such that the line towards zero 
temperature resumes again its high temperature normal 
conducting state linewidth. 
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The shape of the Cu2 lines in both oxygen exchanged 
samples are identical and always asymmetric with a tail 
towards lower frequencies. To improve the reproducibil- 
ity of the analysis we decided to take the center of gravity 
of the upper half of the line as the position of the line. 
The scatter of the line's position as defined above lies 
in the range of ± 2 kHz. Into the same range falls also 
the scatter of the linewidth defined as the full width at 
the half of the line's maximum (FWHM). Figs. and | 
exhibit the observed temperature dependence of the po- 
sition and the linewidth of the Cu2 line in 16 O and 18 O 
YBa2Cu40g samples. At first one clearly notices that 
the frequency data points of the 18 O sample lie beneath 
the ones from the 16 O sample and that the frequency 
shift of about 15 kHz between the two sets of data is 
roughly constant. A closer inspection, however, reveals 
that there is also a slight temperature shift between the 
two temperature dependences. This can be seen best at 
the temperatures where the frequencies have their min- 
ima and the temperature dependences their largest slope. 
These temperatures lie lower for the 18 sample. Fur- 
ther one notices the extraordinarily smooth passage of 
the frequency from the normal to superconducting state. 
There is no jump in the frequency or even a slight change 
in its slope at T c . Also unusual is the rather abrupt sat- 
uration of the temperature dependence of the frequency 
below 30 K. In contrast the linewidth of both 16 and 
18 samples (Fig. ||) show within error no difference in 
the temperature dependence. The dependence itself is 
surprising and has not been observed in high-T c cuprates 
previously. Remarkable is the dramatic change of the 
course of the normal state linewidth temperature depen- 



FIG. 1: Temperature dependence of the NQR frequency of 
plane 63 Cu2 in 16 O (full triangles) and 18 O (full circles) ex- 
changed YBa2Cu4C>8. The solid lines are fits of Eq. (^) to the 
experimental data and the broken lines represent the conjec- 
tured temperature dependence of ^q( 63 Cu2) that comes from 
the thermal expansion of the lattice. 



IV. DISCUSSION 

A. NQR frequency 

With help of vq one is able to monitor the electric 
field gradient (EFG) at the nuclear site of interest. The 
EFG, a ground state property of a solid, depends sen- 
sitively on the charge distribution in the solid. Thus 
knowledge of the EFG and its temperature dependence 
in a high-T c superconductor can contribute valuable in- 
formation concerning the electronic properties of the ma- 
terial. The two naturally occurring copper isotopes 63 Cu 
and 65 Cu have nuclear spin 3/2 and quadrupole moments 
63 Q and 65 Q. In a EFG the Cu nuclei experienc e NQR at 
a single frequency j^ 3 ' 65 = ( 63 ' 65 QeV zz /2h) • yT + ?f/3 
where V zz is the major principal axes of the EFG ten- 
sor and r\ the EFG asymmetry parameter defined as 
V = {V xx — Vyy)/V zz . Here we used the convention: 
\V XX \ < \V yy \ < \V ZZ \. To discuss the results a compar- 
ison with EFG calculations is necessary. The ab initio 
calculations of the EFG |(| are unfortunately not yet 
precise enough for a description of small changes of vq . 
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FIG. 2: Temperature dependence of the linewidth of the plane 
63 Cu2 NQR line in 16 and 18 exchanged YBa 2 Cu 4 08. The 
to the linewidth data in the normal 



solid line is a fit of Eq. 
conducting state. 



We use therefore a semi-empirical approach to the prob- 
lem. We separate the EFG into two contributions: (i) 
the lattice part that refers to the charge distribution of 
the lattice ions surrounding the atom containing the nu- 
cleus in question and (ii) the valence contribution coming 
from electrons (holes) of the incompletely filled electronic 
shells of the same atom. The principal components of 
the EFG tensor V aa add as: V aa = V£g tice + V™ lence , 
where V^ tice and V™ lence denote the lattice and the 
valence part, respectively. Due to the point symmetry at 
the plane-copper (Cu2) site the principal axes of the EFG 
tensor at this site lie along the crystallographic axes with 
the major principal axis V zz parallel to the crystal c-axis 
||. The Cu2 atoms in the Cu02 plane act closely like 
Cu 2+ ions with a hole in the x 2 — y 2 orbit. Bleaney and 
coworkers Q show from electron spin resonance studies 
that a single hole in the d-shell gives rise to an axial EFG 
at the Cu nucleus equivalent to about 70 MHz. The sign 
of Vaa lence from the hole in the x 2 — y 2 orbit is oppo- 
site to the sign of V^ Uce produced mainly by the four 
negatively charged neighbour oxygen atoms. A rough es- 
timate with the so-called point charge model, where the 
lattice ions are approximated by point charges, delivers 
yiattice half the size of V™ lence . Therefore in case of Cu2 
the sum of the two EFG parts of opposite sign turns out 
to be a difference governed by the larger valence part. 
The valence contribution to the EFG remains tempera- 
ture independent (the valence does not change) whereas 
the lattice part diminishes with increasing temperature 
due to the thermal expansion of the lattice. Applied to 



Cu2 this would mean an increase of the summed EFG and 
thus an increase of vq with increasing temperature. Such 
a behavior is indeed observed for Cu2 in YBa2Cu40s, 
but only for temperatures above 200 K. Below 200 K 
the behaviour changes completely so that towards lower 
temperatures vq instead decreasing starts to grow and 
finally saturates at 30 K. Since a high resolution neu- 
tron powder diffraction study |J shows that the crystal 
axes of YBa2Cti408 increase smoothly with temperature 
a change of the vq temperature behavior at 200 K signals 
something new, a feature that most likely is connected 
to the electronic system of the Cu02 planes. That the 
lattice of YBa2Cu40s indeed behaves in a conventional 
way one can also learn from the temperature dependence 
of ^Q(Ba), which smoothly decreases with temperature 
from the lowest (5 K) to the highest (400 K) temperatures 
[jl0|| . The Ba ion has a filled electron shell and therefore 
experiences only the lattice part of the EFG which due 
to the thermal lattice expansion decreases with temper- 
ature. The authors of Ref. pi could fit the tempera- 
ture dependence of vq(Q&) quite accurately by a power 
law: vq{T) = vq(Q)(1 - AT 143 ). In an effort to bet- 
ter discern the effect of the new electronic feature at low 
temperatures we attempt to extrapolate the normal high- 
temperature behaviour of z/q(Cu2) to temperatures be- 
low 200 K. In doing so we assume that the power law 
used for VQ(Ba) holds also for the lattice part of Cu2 
EFG. After subtracting this extrapolation from the orig- 
inal j^q(Cu2) data we get a "rest" whose temperature de- 
pendence looks very familiar, namely like the Cu2 Knight 
shift temperature dependence turned on head jll]] . Lack- 
ing more insight we decide to fit this "rest" by an empir- 
ical function _Btanh 2 (^), that is very reminiscent of the 
empirical spin-pseudogap function used to fit Cu2 Knight 
shift and spin-lattice relaxation data |l], |llj]. In the ac- 
tual ^q(Cu2) vs. T fit we combine the above two steps 
and use the expression: 



v Q (T) = C + AT 



1.43 



Btanh 2 ( — ) 
\ 2T ) 



(1) 



with the constants A, B, C and A as free fit parameters. 
The plus sign of the second term in the expression for 
vq(T) takes care of the fact that in case of Cu2, in con- 
trast to Ba, the smaller lattice contribution to the EFG 
gets subtracted from the larger valence one. The ob- 
tained fit is quite good (see Fig. |l|). The corresponding 
fit parameters are collected in the Table |. The differ- 
ence in the quantum mechanical zero-point displacement 
and thermal lattice expansion can explain the nearly con- 



TABLE I: The parameters A, A, B, and C used in EqjT]) to 
fit the temperature dependent NQR frequency of plane Cu2 
in 16 and 18 exchanged YBa 2 Cu 4 8 . 



A (K) 


A (MHzK- 1 ' 43 ) 


B (MHz) 


C (MHz) 


lb O 188.0(1.6) 
ls O 180.0(1.6) 


2.59(7) xl0" b 
2.69(6) xl0~ 5 


0.184(2) 
0.182(2) 


29.642(3) 
29.629(2) 
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stant 15 kHz shift in the lattice part of the EFG be- 
tween the 16 and 18 samples. The heavier 18 has 
a smaller zero-point and thermal fluctuation amplitude 
than 16 and therefore expands the lattice less than 16 
as observed by the x-ray measurement of the lattice pa- 
rameters. Consequently the lattice part of the EFG in 
the 18 O sample is larger than in the 16 O sample which 
leads in case of Cu2 to a lower vq. Of greater interest 
is the parameter A representing the energy scale of the 
vq feature. Its value though close to the magnitude of 
the spin-pseudogap is definitely smaller than the later 
one. Further, the oxygen isotope effect on A is substan- 
tial. The corresponding partial oxygen isotope effect co- 
efficient is oivq — 0.42(11) and thus much larger than 
both the spin-pseudogap coefficient ape = 0.061(8) and 
the T c coefficient a Tc = 0.056(12) [0. The EFG does 
not depend on spin but exclusively on charge, however, 
on all charges, irrespective what state they occupy. In 
that respect the EFG differs from the Knight shift and 
spin-lattice relaxation which experience only the spins of 
the charge carriers whose energies are close to the Fermi 
level. Therefore, in case of the EFG, at the formation 
of a pseudogap all the electrons have to be considered, 
not just those close to Fermi surface. This makes conclu- 
sions on the basis of an empirical analysis as ours rather 
uncertain. Nevertheless, the observed large oxygen iso- 
tope effect on the temperature dependence of ^q(Cu2) 
indicates that the involved charge feature in the CuC-2 
planes is influenced by the coupling of the charge carri- 
ers to the lattice as it is the case for spin-pseudogap and 
superconductivity. 



B. Linewidth 

The linewidths of the Cu2 NQR lines from the two 
oxygen isotope samples (Fig. |^) are equal within error 
at all temperatures. Further, we find that at 350 K 
the linewidth ratio of the 63 Cu and 65 Cu isotope NQR 
lines is 1.085(7) which within error equals to the ratio 
(1.0806) of the isotopes' nuclear quadrupole moments. 
This allows the conclusion that at 350 K the line is broad- 
ened predominantly by quadrupolar effects produced by 
static EFG inhomogcncitics in the material. Such in- 
homogeneities, in principle, could be of intrinsic ori- 
gin generated for instance by a charge instability in the 
electronic system of the plane. However, this possibil- 
ity most probably can be excluded since a comparison 
of linewidths measured at 350 K in YBa2Cu40s sam- 
ples from different batches, shows a variation of these 
quadrupolar linewidths depending on the parameters of 
material preparation. Since the superconducting prop- 
erties are not sensitive to the structural disorder there 
is not much incentive to invest for NQR purpose alone 
into the time consuming improvement of the prepara- 
tion parameters necessary to make a structurally perfect 
material. Nevertheless, the obtained 63 Cu2 line with its 
125 kHz linewidth at room temperature is to our knowl- 



edge the narrowest plane 63 Cu NQR line yet observed in 
high-T c cuprates. Though this proves the high quality 
of our samples improvements certainly are still possible, 
especially for a stoichiometric material as YBa2Cu40s 
where one expects an order of magnitude narrower lines. 
The observed lines are broadened asymmetric having a 
tail towards lower frequencies. Comparison with the Ba 
NQR line suggests that this asymmetry has to come from 
a lattice imperfection that produces simultaneously an 
EFG inhomogeneity at Cu2 and Ba sites since the shape 
of Ba NQR line is asymmetric too but with a tail to- 
wards higher frequencies. Lattice defects can produce 
various EFG inhomogcneities. The defect we observe, 
is somewhat special since it can produce opposite line 
asymmetries for Cu2 and Ba NQR lines. A possible can- 
didate that could do so would be for instance a lattice 
defect that locally shrinks the lattice and thus increases 
the lattice part of EFG at both Cu2 and Ba sites. As 
already mentioned, any increase of lattice EFG pushes 
Ba and Cu2 NQR frequencies in opposite direction, it in- 
creases the Ba and decreases the Cu2 one what explains 
the opposite lineshape asymmetries. The asymmetry of 
the lines does not change with temperature. Even though 
asymmetric, the Cu2 line is narrow enough to allow a 
very precise measurement of the broadening of the line 
with decreasing temperature in the normal conducting 
state and what is more important it enables in a high T c - 
cuprate to discern how a plane-Cu line gets narrow again 
below T c . The observed change in the temperature be- 
havior at T c is extremely sharp with a very rapid decrease 
of the linewidth at the passage into the superconduct- 
ing state where the line continues to narrow to approach 
at zero temperature its smallest normal state high tem- 
perature value. This behaviour, most likely intrinsic in 
YBa 2 Cu40s, has not been observed yet in other high-T c 
cuprates. It is the excessively broad plane-Cu lines in 
other cuprates that do not allow a similar observation. 
Most of cuprates are plagued by structural disorder gen- 
erated by the only possible nonstoichiometric doping al- 
ternative. In YBa2Cu40g at temperatures just above T c 
where the Cu2 line is broadest the ratio of the linewidths 
of the 63 Cu and 65 Cu isotopes decreases from the high 
temperature value to 1.06(1) which is smaller than the ra- 
tio of the two isotopes' quadrupole moment and thus sig- 
nals very likely a magnetic component in the line broad- 
ening. In case of pure magnetic broadening of the line 
the ratio of the linewidths of the two isotopes would be 
equal to the quotient of the corresponding Cu isotope 
gyromagnetic ratios e3 j/ 65 j = 0.9335. That a magnetic 
broadening of the line with decreasing temperature in- 
deed takes place we can also infer from the temperature 
dependence of the indirect Gaussian contribution 1 /T^q 
to the Cu2 spin-spin relaxation rate (T^| that comes from 
the nuclear magnetization transfer via the itinerant elec- 
tron spin system and which is proportional to the static 
electron spin susceptibility (x(Qaf,u = 0)) at the anti- 
ferromagnetic wave vector Qaf- With decreasing tem- 
perature x{Qaf,u = 0) exhibits a Curie like behaviour 
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Jl2{ . Its contribution to spin-spin relaxation is Curie like 
and presumably the same is true for the linewidth. To 
find out the temperature dependent contribution to the 
linewidth we try to deconvolute the total linewidth into 
its temperature independent and dependent parts. We 
find the linewidth at the lowest temperatures to be equal 
to the linewidth at the highest temperature where the 
magnetic contribution due to the Curie like behavior of 
x{Qaf,u — 0) is expected to be rather small and the 
only remaining contribution is then quadrupolar as we 
already noticed before. Usually a quadrupolar contri- 
bution that comes from lattice imperfections and has the 
observed large size does not vary much with temperature. 
Since we find the linewidths at the two temperature ex- 
trema equal we assume that we have an "underground" 
quadrupolar contribution that remains constant through 
out the whole temperature range. In general the broad- 
ening of a line caused by a combined operation of two 
distinct mechanisms is rather complex. There are few 
cases where simple relations exist as for instance for two 
Lorentzian or two Gaussian like line broadenings where in 
the first case the contributions add linearly and in the sec- 
ond quadratically. We do not expect that the broadening 
we study has one of these simple forms. However, to be 
able to proceed we make an assumption that the constant 
and the temperature dependent part of the linewidth add 
with equal power n. We keep n as a free fit parameter 
with a value that lies somewhere between 1 and 2. We 
use the linewidth data from the normal conducting state 
as a gauge to find out the temperature dependent part 
of the linewidth in the superconducting state. For this 
purpose we fit the data gained in the normal conducting 
state so that we keep the underground quadrupolar part 
of the linewidth constant and allow the temperature de- 
pendent part to vary Curie like with temperature. For 
the fit procedure we use the following expression: 



(6v t ot(T)) n = {8v quad ) n + (C/T) n . 



(2) 



The fit of the normal conducting state data yields: 
5v q uad = 124(1) kHz, C = 3.0(4) MHzK, and power 
n = 1.24(7). With help of the parameters 8v qua d and 
n we afterwards decompose the measured linewidth data 
to get the wanted temperature dependent contribution 
to the linewidth, 8v(T). The final result is presented 
m Fig. § As one can see the temperature dependent 
part of the linewidth tends to a value close to zero when 
the temperature approaches zero. The temperature de- 
pendent part of the linewidth can be of magnetic and 
of quadrupolar origin. Due to the large static electronic 
susceptibility at Qaf the magnetic part of the linewidth 
comes predominantly from the staggered magnetization 
induced by intrinsic as well as extrinsic magnetic field 
components at Qaf- From l/T^f measurements in the 
superconducting state Q we know that x{Qaf,^ = 0) 
decreases only 15% of its value at T c when going far into 
the superconducting state. Therefore, the reduction of 
the magnetic part of the NQR linewidth in the supercon- 
ducting state has to come from a reduction or screening 



of the extrinsic magnetic fields in the superconducting 
state. From our measurements it is obvious that with 
decreasing temperature the NQR line broadens smoothly 
as long as the sample stays normal conducting but starts 
sharply to narrow when the sample turns superconduct- 
ing. The whole increase in linewidth accumulated from 
350 K down to T c disappears away in the superconduct- 
ing state. At the moment we do not have an adequate 
explanation of this very unusual behavior of the linewidth 
of the Cu2 NQR line below T c . 
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V. SUMMARY 

We performed accurate measurements of the temper- 
ature dependence of the 63 Cu2 NQR line frequency and 
linewidth in normal and superconducting 16 and 18 
exchanged YBa2Cu 4 0g. At the transition of YBa2Cu40g 
from normal into the superconducting phase we observe 
that ^q(Cu2) passes T c very smoothly without a disconti- 
nuity either in value nor in slope. In contrast to ^q(Cu2), 
the linewidth exhibits a drastic change of the course at T c . 
While increasing in the normal conducting phase down to 
T c , the linewidth starts to decrease sharply in the super- 
conducting phase to finally resume towards zero temper- 
ature its high-temperature value in the normal conduct- 
ing phase . Such a behaviour of the narrow plane-copper 
NQR line in YBa2Cu40s has not been observed yet in 
other high-T c cuprates. At the moment we do not have 
an adequate explanation of this extraordinary tempera- 
ture dependence of the Cu2 linewidth. The frequency 
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vq (Cu2) of 18 exchanged YBa2Cu40s is always lower 
than vq (Cu2) of 16 exchanged YBa 2 Cu40 8 . The nearly 
constant frequency shift of 15 kHz can be accounted for 
by the difference in the quantum mechanical zero-point 
displacement and thermal lattice expansion between the 
16 and 18 samples. More important, there is a well dis- 
cernible temperature shift between the ^q(Cu2) temper- 
ature dependencies of the 16 and 18 samples. Lacking 
a detailed theoretical description we attempt empirically 
to decompose the temperature dependence of ^q(Cu2) 
into two parts: one part coming from the thermal ex- 
pansion of the lattice and the other produced by the 
charge redistribution during the formation of new elec- 
tronic structures in the CuC>2 planes. From the fit of the 
empirical expression for vq(T) to the experimental data 



we determine for the conjectured formation of new elec- 
tronic structures an energy scale A = 188.0(1.6) K for the 
16 and 180.0(1.6) K for the 18 exchanged YBa 2 Cu 4 8 . 
The corresponding partial oxygen isotope effect coeffi- 
cient a UQ = 0.42(11) is much larger than both the spin- 
pseudogap coefficient ape — 0.061(8) and the T c coeffi- 
cient u Tr = 0.056(12) §. 
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